Thermal conductivity of Mg-doped CuGeOs at very low temperatures: Heat 
conduction by antiferromagnetic magnons 
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Thermal conductivity k is measured at very low temperatures down to 0.28 K for pure and Mg- 
doped CuGeC>3 single crystals. The doped samples carry larger amount of heat than the pure sample 
at the lowest temperature. This is because antiferromagnetic magnons appear in the doped samples 
and are responsible for the additional heat conductivity, while n of the pure sample represents 
phonon conductivity at such low temperatures. The maximum energy of the magnon is estimated 
to be much lower than the spin-Peierls-gap energy. The result presents the first example that k at 
very low temperatures probes the magnon transport in disorder- induced antiferromagnetic phase of 
spin-gap systems. 

PACS numbers: 66.70.+f, 75.30. Kz, 75.50.Ee 



Recently, impurity-substitution effect on the spin- 
singlet ground state has been intensively studied in a 
variety of low-dimensional spin systems, and the results 
indicate that only a slight substitution of non-magnetic 
impurity essentially changes the ground state. When 
the disorder is introduced by the non-magnetic impu- 
rity, antiferromagnetic (AF) ordering immediately ap- 
pears without destroying the spin-gap feature, .commonly 
in a spin-Peierls (SB} compound CuGe03,ETEl a ladder 
compound SrCuaQ3,B and a newly found Haldane sys- 
tem PbNi2V2O8.0tl Since it was believed that spin-singlet 
state and antiferromagnetism are mutually exclusive, the 
case of CuGeOs, in which both the remaining SP or- 
dering and the AF ordering developing -with impurity- 
doping turned out to be long-range order ,0 is remarkable. 
Fukuyama et al. proposed a theoretical model for this 
anomalous j)hase, i.e., dimerized antiferromagnetic (D- 
AF) phaseB In their model, staggered AF moments and 
lattice dimerization are spatially modulated, keeping the 
AF and SP long-range correlations; degree of the lattice 
dimerization is larger in between the impurity sites than 
around the impurity sites and magnitude of the staggered 
moments are larger near the impurity sites. 

Even though the coexistence of the two long-range or- 
ders are established in lightly-doped CuGeC>3, it is to be 
elucidated how low-energy spin excitations of SP state, 
which holds spin gap, and those of AF state, which is 
gapless in isotropic systems, compromise with each other 
in D-AF phase. Recently, Saito and Fukuyama extended 
the theory of Ref. [)] and predicted that a "slow" gap- 
less AF-magnon branch shows up in addition to the 
gapped SP mode when impurity is doped.L3 They pre- 
dicted that the energy scale of the magnon branch is 
smaller than the SP gap (~ 24 K). Existence of this in- 
gap magnon branch is observed by ueuttpn scattering 
for rather highly substituted samples JlJil 3 ] (Zn concen- 
tration was 3.2% for the sample used in Ref. [fl], for ex- 
ample). However, since recent study using high-quality 
Mg-doped single crystals, in which impurities distribute 



more homogeneously than in Zn-doped systems, has re- 
vealed disorder-induced first-orcles transition around the 
Mg concentration of x c — 2.4%,EJ onbz-below which the 
long-range SP ordering is established,!!^ it is essential to 
examine the magnetic excitations of samples with impu- 
rity concentration less than x c . In this work, we have 
measured thermal conductivity k at very low temper- 
atures, using lightly Mg-doped samples. The samples 
are cooled down to 3 He temperatures, in order to study 
the low-lying excitations. The thermal conductivity of 
the doped samples exceeds that of the pure sample at 
the lowest temperature. We will show that the in-gap 
magnons, which is intrinsic to D-AF phase, indeed exist 
and are responsible for the excess low-temperature heat 
transport in the doped samples. 

The Cui_ x Mg x Ge03 single crystals were grown with 
a floating-zone method. The Mg-concentration is de- 
termined by inductively coupled, plasma-atomic emis- 
sion spectroscopy (ICP-AES).EjO For the thermal con- 
ductivity measurement, we use pure, x = 0.016, and 
x = 0.0216 samples, all of which were already well char- 
acterized using dc susceptibility and synchrotron x-ray 
diffraction measurements.Efrtlj The transition tempera- 
tures are shown in Table [j] for each sample. SP long- 
range order is observed at T$p as a resolution limited 
FWHM of the x-ray Bragg peak from lattice dimeriza- 
tion. The Neel tempeeature TV is determined by the 
magnetic susceptibility.tJ 

Thermal conductivity is measured down to 0.28 K with 
3 He refrigerator using "one heater, two thermometers" 
technique. Gold wires, which are tightly connected with 

TABLE I. 



Sample 


X 


T$» [K] T N [K] 


w [mm] 


<n [ m /s] 


A 





14.5 


0.17 




B 


0.016 


10.5 2.5 


0.17 


70 


C 


0.016 


10.5 2.5 


0.24 


70 


D 


0.0216 


8.5 3 


0.10 


140 



1 



^ 1.5 
E 

S 1 

* 0.5 


0.25 



III 

: C Ul _ Mg Ge0 3 


iii 

Sample A / 
(x = 0) : J 

N. « _ 
\ # 


Sample B 
(x = 0.016) 




, , ^•■r , , , 


^ (a)" 



0.2 0.4 0.6 0.8 1 1.2 

T[K] 



^ 0.2 
E 0.15 

i o.i 

0.05 




Sample B 

(x = 0.016) y 




Sample A 
(x = 0) 



(b) 







0.05 



0.15 



0.2 



0.1 
T 3 [K 3 ] 

FIG. 1. (a) Thermal conductivity of Cui-^Mg^GeOa sin- 
gle crystals below 1.3 K. Solid curve is for one of the x = 0.016 
samples (B) and dashed curve is for the pure sample (A), (b) 
Thermal conductivity of the same sample as a function of T 3 . 
The dashed lines show that n is proportional to T [i below 
0.58 K for Sample A. 



a microchip heater and two calibrated RuO sensors, are 
attached on the samples by GE-varnish. Temperature 
difference between the two thermometers are typically 
3% of the sample temperature. Since we will discuss 
the thermal conductivity in the low-temperature limit 
(Casimir's limit), EZHl3 where heat carriers are scattered 
dominantly by crystalline boundaries, we paid special at- 
tention both to the sample size and to the smoothness of 
the boundaries. Typical dimension along the c-axis is 
around 3 mm and thermal gradient is applied along the c 
axis. As shown in Table |, the geometrical mean widths w 
(square root of the cross section) , which is proportional 
to the mean free path in the Casimir's limitjL3 of the 
pure sample (Sample A) and one of the x — 0.016 sam- 
ples (Sample B) are set identical for direct comparison in 
k. The smooth boundaries of the crystal is achieved by 
cleaving for the b-c surfaces and by cutting with a sharp 
razor blade for the a-c surfaces. Also, we have measured 
specific heat C, in order to derive the mobility of the 
relevant heat carriers. The specific heat measurement is 
carried out down to 0.4 K with the commercial PPMS 
(Quantum Design) heat-capacity probe, using a relax- 
ation method. The mass of the samples is typically 3 mg. 
The samples used for the specific heat and the thermal 
conductivity measurements are cut from the same piece 



of the crystals. 

Figure 1(a) shows temperature dependence of k for 
Samples A and B below 1.3 K. In Fig. 1(b), k in the tem- 
perature range below 0.58 K (T 3 < 0.2 K 3 ) is plotted as 
a function of T 3 . k of the pure sample rapidly decreases 
with decreasing temperature and becomes proportional 
to T 3 below 0.58 K. At 1.3 K, k of the x = 0.016 sample 
is smaller than that of the pure sample, so that our prevL. 
ous results in higher temperature range are reproduced.EEl 
It is natural that k is suppressed in the presence of im- 
purities owing to the scattering by disorders. However, k 
of the x — 0.016 sample exceeds that of the pure sample 
below 1 K down to the lowest temperature, which is not 
explained by the above simple picture. 

First, let us discuss k of the pure sample. Magnetic ex- 
citations are negligible in the pure sample below 1.3 K, 
because the spin-gap energy of CuGeC>3 is more than, 
one order of magnitude larger than the temperatures. EHI 
Therefore, the heat is dominantly carried by phonons 
there. Assuming kinetic approximation, the phonon ther- 
mal conductivity K p h is written as 



K ph — 77 ^-"ph^ ph^ph, 



(1) 



where C p h is specific heat, i> p h is velocity and l p h is mean 
free path of the phonons. Since C p h also depends on 
temperature as T 3 at such low temperatures, lph should 
be independent of temperature below 0.58 K. This result 
means that the phonon conductivity reaches the Casimir 
limit, where the mean free pathis. determined simply by 
the dimension of the crystal £ZEaE!J For a rectangular- 
shaped crystal l p h is given as 

l ph = 1.12<£, (2) 

assumimg isotropic phonons.EEl We have measured the 
specific heat independently and examined the validity 
of Eqs. (1) and (2). First, u p h is calculated from k 
and C data by EqSj— ,(1) and (2), as v p h = (3k) /(C ■ 
1.12u)) ~ 1600 m/s.E 3 ] On the other hand, t> p h can be 
estimated only from the low-temperature C data, as- 
suming the Debye model. Thus obtained value of v p h is 
~ 1800 m/s, which is in good agreement with the estima- 
tion from the k data. The result shows that the relation 
of Eqs. (1) and (2) is satisfied at the temperatures shown 
in Fig. 1(b) for the crystal used. 

Noting that phonon conductivity governed by the 
boundary scattering gives the maximum value of K p h (any 
additional scattering would suppress the conductivity), 
we can notice that K p h of Sample B cannot exceed n of 
Sample A in Fig. 1(b), using Eqs. (1) and (2), because 
w of Sample A and B is identical (Table |) and little 
x dependence is expected for C p h and i; p h- Therefore, 
the result of the excess k in Sample B requires an addi- 
tional excitations which can carry heat at low tempera- 
tures down to 0.28 K in the Mg-doped sample. Consider- 
ing the AF ordering in the impurity-doped CuGeOs, it is 
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most likely that antiferromagnetic magnons are responsi- 
ble for this excess low-temperature heat conductivity in 
the Mg-doped sample. 
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FIG. 2. Specific heat and thermal conductivity of 
Cui-zMg^GeOs with the concentration of x = 0.016 below 
~ 0.58 K as a function of T 3 . Specific heat of the pure 
CuGcC>3 is plotted together. Inset: k/w of two x = 0.016 
samples with different w, i.e., Samples B and C. 

In order to examine whether the Mg-doped sample 
also satisfies the Casimir's condition, under which more 
quantitative discussion is possible, we compare low- 
temperature k of Sample B and that of another x = 0.016 
sample with different w (Sample C). The upper inset of 
Fig. 2 shows that the values of K are different between 
Samples B and C. On the other hand, as shown in the 
lower inset, the low-temperature k/w data (T 3 < 0.1 K) 
of the two samples show only 3% difference, which is 
mainly due to the error in determining the distance be- 
tween the gold wires connected to the thermometers. 
Since the result means that k differs in proportion to 
w at low temperatures, it is strongly suggested that the 
Casimir's condition is satisfied for both samples and that 
both K p h and the magnon heat transport governed 
by the boundary scattering there, i.e., 



= ^C m v m 1.12a c w 



(3) 



[a c (= \Jv m v m /v m 2 ~ \/\JaJb\/Jc) ijHa factor due to 
the anisotropy of the magnon velocityta and v m is the 
magnon velocity in the i-th-axis direction). 

Assuming Eq. (3), we can crudely estimate v m with the 
value of C m obtained by the independent specific-heat 
measurement and examine whether the magnon branch 
is in the SP gap. The specific heat of Samples A and B 
is shown in Fig. 2 as a function of T 3 in the tempera- 
ture range below - 0.58 K 3 (T 3 < 0.2 K). Since C of 
the pure sample (diamonds in Fig. 2), which represents 
C p h in this temperature region, is more than one order 



of magnitude smaller than that of the Mg-doped sam- 
ples, we can neglect C p h and assume C rn ~ C for the 
doped sample. K m of the x = 0.016 sample is estimated 
by subtracting k of the pure sample and is plotted to- 
gether in Fig. 2. The K m (T) data show exactly the same 
temperature dependence as the C(T) data below 0.46 K 
(T 3 < 0.1 K), indicating again that Eq. (3) is satisfied in 
this temperatupj range. Assuming J c ~ 10 Jb ~ 100| J a \ 
and a c ~ 0.1, EJ v m is estimated approximately to be 
70 m/s from Eq. (3). The AF-magnon energy is lower 
than q c v m ~ 0.12 meV [q c (= 7r/4c) is the wave number 
at the center of the Brillouin zone and c (= 2.9 A) is 
the distance between adjacent Cu atoms along the Cu-0 
chain] . Since this value is two orders of magnitude lower 
than the SP gap (A$p ~ 2 meV), we can conclude that 
the additional heat transport observed in the Mg-doped 
CuGeC>3 is due to the in-gap AF-magnons.Eil 
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3. (a) Thermal conductivity of Cui-zMg^GeOa 
(x — 0, 0.0216) crystals, divided by w is plotted against T 3 . 
(b) Specific heat and thermal conductivity of Cui^^Mg^GeOa 
(x = 0.0216) below ~ 0.2 K as a function of T 3 . Specific heat 
of the pure CuGeOs is plotted together. 



The magnon velocity can be estimated also for the 
x = 0.0216 sample (Sample D) in the same manner. In 
order to subtract K p h, it is convenient to compare k/w 
of Sample D to that of Sample A, which represents the 
phonon contribution, k/w of Samples A and D is plot- 
ted against T 3 in Fig. 3(a). k/w of Sample D is larger 
than that of the pure sample. The difference corresponds 
to the magnon contribution K m /w of Sample D. In Fig. 
3(b), K m /w and independently obtained C data are plot- 
ted together. Following the same discussion as that for 
the x = 0.016 sample, v m is estimated approximately 
to be 140 m/s. The upper limit of the magnon energy 
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q c Vm is ~ 0.23 meV, which is twice as high as that of the 
x = 0.016 sample but is still much lower than the Agp. 

The velocity of AF-magnon for the usual uniform Neel 
state is given by 2zJ c Sc^B where z is a factor of or- 
der unity, J c is the interaction energy and S is the 
spin in each magnetic site. Even when quantum fluc- 
tuations, which, diminish the effective value of S to 0.2 
times smaller,LJ is taken into account, the magnon ve- 
locity of as fast as ~ 1000 m/s is expected, (assuming 
J c = 120 K).cJ Since we obtained much smaller v^, it 
is shown that effective S is significantly suppressed ow- 
ing to the SP ordering. In the model of Refs. ^ and 
|l0| , the staggered moment is strongly suppressed in be- 
tween the impurity sites. As the result, a sort of spatially 
"averaged" spin, which is directly proportional to the 
magnon velocity, becomes much smaller than 1/2. Our 
crude estimation from the low-temperature k gives v c m of 
the x = 0.0216 sample approximately twice as large as 
that of the x = 0.016 sample, indicating that rapidly 
increases with x. Such x dependence is consistent with 
the calculation in Ref. [l0| Note that of around 1300 
m/s can be estimated for the Zn-3.2%-doped sample from 
the neutron datajiil which is nearly one order of magni- 
tude larger than the value estimated for our x = 0.0216 
sample. 

It should be emphasized that the mean free path 
of the AF-magnon in the SP gap reaches as long as 
I.12a c w ~ 18 /im, according to Eq. (3), i.e., the magnons 
are mobile to a distance 60,000 times longer than the 
spin-distance (c) without being scattered. Such coherent 
motion of the magnons is only possible with extremely 
long correlation length of the corresponding magnetic 
order SB Therefore, the mixture of SP and AF ordering in 
the doped CuGeC>3 is certainly a true long-range order. 

In summary, we have measured thermal conductivity 
of pure and Mg-doped CuGeC>3 single crystals at very 
low temperatures, in order to examine the anomalous 
low-temperature phase in the impurity-doped CuGe03, 
where SP and AF order coexists. While the low- 
temperature thermal conductivity is dominated by K p h 
in the pure CuGeOs, magnons also carry considerable 
amount of heat in the Mg-doped samples. Estimating the 
magnon velocity, we have shown that the AF-magnons 
are present in the SP gap, as predicted by Saito and 
Fukuyama. It is demonstrated that thermal conductivity 
is a powerful tool in elucidating the low-energy magnetic 
excitations in the disorder-induced AF phase of spin-gap 
systems. Our next step is to examine the excitations in 
the uniform AF phase in highly substituted CuGeOs, and 
to seek differences from that of D-AF phase. 
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in- Aid for COE Research of the Ministry of Education, 
Science, Sports, and Culture. 
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